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ABSTRACT 



Vertical stern accelerations for five geometri- 
cally similar destroyer hull forms of different lengths 
in the 250-550 foot range are predicted on the basis of 
model tests for a ship speed of 17 knots at various 
headings in representative short-crested irregular seas 
of varying severity. 

Using the concept of statistical inference, the 
model test predictions are treated as a statistical 
sample and used to predict what would occur in the total 
population of possible stern accelerations. 

Then, with strong emphasis on the systems 
approach, it is shown how these predictions can be used 
as input to a probability analysis that leads to a pre- 
diction of long range system behavior/effectiveness, 
the system under consideration being the Destroyer/ 
Variable Depth Sonar System. 
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N OIvTEN C LAT1JHE 



DD - Destroyer type ship. 

VDS - Variable Depth Sonar. 



RAO - Resoonse Amplitude Operator. 

- Actual wave length of regular waves in 
feet measured perpendicular to lines of 
crests . 

I 

Effective - Wave length apparent to a ship proceeding 
obliquely through regular waves. 

L - Ship length in feet measured along 

designer's waterline 



Amplitude 



He ight 



X 



So 






- For a regular or irregular function, 
peak to peak value civided by two. 
Sometimes referred to as single 
amplitude . 

- For a regular or irregular function, 
peak to peak value. Sometimes referred 
to as double amplitude. 

- Amplitude of ship pitch angle in radians; 
vertex of angle is defined as being at 
ship's center of gravity. 

- Instantaneous value of pitch. 

- Amplitude of heave (vertical transla- 
tion) of ship in feet at center of 
gravity. 

- Instantaneous value of heave. 



0 - Amidships 

2 X - Vertical translation at a point x feet 

from 0. 

h - Wave amplitude in feet. 
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h 


- Instantaneous value of distance in feet 
between wave surface and calm water zero 
datum* 


\}"max 


- Maximum wave slope. See Page 9. 


K 


- Phase angle between pitch and wave* 
See Pages 8 and 9. 




- Phase angle between heave and wave. 
See Pages 8 and 9* 


Kr 


- Phase angle between heave and pitch* 
See Page 9. 


^ /%, A-> 


- See Figure 3 on page 11 for angular 
nomenclature * 


V 8 


- Ship speed in feet per second. 


to 


- Circular frequency of wave in radians per 
second . 




- Encounter frequency in radians per second. 


t 


- Time in seconds. 


g 


- Acceleration due to gravity; 32.2 ft. per 
(second )* 


X 


- Distance along ship's length from 8 to 
some point, (+) forward, (-) aft. In 
Appendix A, ordinate of irregular stern 
acceleration record. 


^X 


- Amplitude of vertical acceleration at a 
point x along ship's length. In this 
thesis x = -181 feet for the 370 foot 
ship length, (x changes with geometric 
similarity for other ship lengths) 


In 


- Logarithm to the natural base e. 


f 


- Increment of. 


corr.f t 2 


- Pierson-Moskowitz sea spectrum ordinate. 
"Corr" implies corrected. See Page 16 


H 


- Lag number; sea spectrum abscissa fre- 
quency coordinate* See Pages 16 and 17, 


A 
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H l/3 


- Significant wave height in feet; defined 
as the average of the l/3 highest wave 
heights extant in any given seaway. 


RSC 


- Response Spectrum Component. See Figure 9. 


IRS 


- Integrated Response Spectrum. See Figure 9. 


1/E 


- Ship response parameter. See Appendix A. 


Sample 


- A relatively small portion of the set of 
values in the universe of a random variable. 


X 


- Mean of sample probability distribution of 
stern vertical "ME acceleration. 


s 2 


- Variance of sample probability distribution 
of stern vertical i(e acceleration. 


3 


- Deviation of sample probability distribu- 
tion of stern vertical Me acceleration. 


N(x,s 2 ) 


- A normal (Gaussian) probability distribu- 
tion function whose mean is x and whose 
variance is S 2 . 


Population 


- The set of all possible values in the 
universe of a random variable. 


/V 


- Mean of population probability distribution 
of 8 tern vertical ME acceleration. 


<r 2 


- Variance of population probability distri- 
bution of stern vertical ME acceleration; 
when used in a general sense , the variance 
of any normal probability distribution 
f unot ion. 


<r 


- Deviation of population probability distri- 
bution of stern vertical Me acceleration. 


V c 

i i 


- True population parameters. 

- Probability of [ ] . 


Au, u 


- Estimate of S'o, ^ o • See Figure 17 on Page 56 


p 


- Location of x and s in parameter space. 


n 


- Size of a statistical sample. 
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N - Number of cycles of a certain ship motion. 

R - Parameter of the Rayleigh probability dis- 

tribution function. 

RMS - Root mean square of some variable 0 See 

Appendix A. 

I/rF - See Appendix A. 



Superscripts 



- Denotes amplitude of a sinusoidal 
function as distinguished from instan- 
taneous value. In the case of an irregular 
function, denotes peak to peak divided by 
two. This superscript does not imply the 
above when used in x; in that case it 
indicates an arithmetic mean. 

- Average of. 



Subscripts 

9 - Measured at ship’s center of gravity. 

max - Maximum value. 

i - heave relative to v/ave . 

y - pitch relative to wave. 

2 V - pitch relative to heave. 

x - at location x feet along ships length 

from B, 

e - Encounter. 

^ - the 4 th value of a variable . 

u - upper. 

L - lower. 

PTP - peak-to-peak value (double amplitude or 

height ) 
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l/Z> - the highest l/3 of a group of values 

of a variable . 

1 - original value, before being changed. 

2 - value after having been changed from 

original value. 
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INTRODUCTION 



The basic method used in this thesis to predict 
the various modes of response of a particular ship to a 
given seaway was introduced over ten years ago (1)* and 
has been extensively employed in the intervening time* 
The challenge of the present day is how to most effec- 
tively utilize such response data in order that the 
systems designer might be provided with reliable long- 
range estimates of system behavior and as well, that the 
system user might have reliable estimates of system 
effectiveness* 

The approach developed in this thesis is a 
broad and general one that is applied to one particular 
aspect of behavior of one system* It is suggested that 
the reader keep in mind the generality of the method, 
while considering the specific case which is developed 
herein for illustration as well as for its current 
naval interest* 

The system considered in this thesis is the 
DD/VDS system. The particular aspect of system behavior 
considered is the vertical acceleration of the YDS tow- 
point at the destroyer *s stern* 

^Numbers in parentheses refer to references listed on 
Page 91* 
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The development undertaken herein can be divided 
into three principal areas, namely: 

(1) prediction of vertical stern acceleration 
(in terms of the 'viB parameter which is 
discussed in detail in Appendix A), for 
various ship lengths, headings and sea 
states ; 

(2) utilization of the concept of statistical 
inference to predict the population dis- 
tribution of vertical stern acceleration* 

(in terms of the *\(E parameter), treating 
the predicted values in (1) above as a 
statistical sample; 

(3) a demonstration of how these predicted stern 
accelerations can be used in a probability- 
based system analysis leading to long- 
range prediction of DD/VDS system behavior 
over the lifetime of such a system. 



*The term " population* is used here to signify all the 
possible values of vertical stern acceleration (in terms 
of the VE parameter) that a destroyer of the particular 
type considered might experience in its lifetime. In 
theoretical probability and statistics, "population” and 
"universe" are used interchangeably to describe the set 
of all the possible values of a given random variable. 



2 



PREDICTION OF STERN ACCELERATION RESPONSE 



General Procedure for Prediction of Ship Response 

In 1953, St. Denis and Pierson (1) introduced a 
method utilizing the principle of superposition whereby 
a ship's response in irregular waves may be predicted 
provided the response in regular waves is known, as it 
is, for example, from model tests in regular waves. 

The method basically consists of (a) representing 
the sea as an energy spectrum, (b) transforming model 
response data taken in regular waves into a response 
amplitude operator (RAO) (referred to by some as a 
transfer function) and (c) multiplying the RAO by the 
sea spectrum to arrive at a response spectrum having 
certain probability properties regarding the response 
being considered. 

One section each will be devoted to the descrip- 
tion of developing RAO's from model data, to the selection 
of sea spectra, ana to the computational procedure used 
to combine sea spectrum and RAO to arrive at a response 
spectrum. 

Development of RAO's from Model Data 

The most comprehensive source of destroyer model 
motions data available at the present time is the Dutch 
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data that can be found in (2). It is particularly 
unique in that model response at various headings to 
the waves was measured and this is, of course, a 
necessary prerequisite if one desires to make predic- 
tions in short-crested seas and to investigate the 
general effect of heading on ship response. This Dutch 
data is also unique in that phase angles were measured 
and presented. The primary limitation of this data is 
that the above variations were done at only one speed, 
namely 17 knots full scale. 

The vertical acceleration at the stern of a ship 
is caused primarily by pitching and somewhat less by 
heaving. Roll is the only other possible contributor 
to vertical acceleration, but its effect was ignored 
because the stern VI'S towpoint is on the ship’s center- 
line. The vertical acceleration also is a function of 
the phase angle between heaving and pitching. Conse- 
quently, the data from (2) used in this investigation 
consisted of heave, pitch and the heave -to-wave and 
pitch-to-wave phase angles. Samples of the model data 
are shown in Figures 1 and 2. 

The model used was 17002. It is a model of 
Her Netherlands Majesty’s Ship Friesland. Its hull 
parameters are as shown in Table I. 
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TABLE I 



L 370 ft. 

Beam 38.6 ft. 

Draft 13.2 ft. 

Blook Coefficient 0.566 



Waterplane Coefficient 0.805 
Priematic Coefficient 0.688 
Displacement 3 S 010 tons 

Number of Rudders 2 



It can be seen from Figures 1 and 2 that a certain 
amount of extrapolation of the model data was necessary. 

On the high end of the ^/L scale, Z Q ->h and ^^211 h/a. = 
ijmax, but the question immediately arose as to Just 
where on the a/L scale it would be reasonable to assume 
that this occurs. Extensive examination of other destroyer 
model data (3), (4) showed that a 'A/L of 3.75 to 4.00 is 
reasonable for heave and pitch having attained their 
limiting values. Since the angle between ships course 
and wave direction had to be taken into account, T 0 was 
faired to effective wave slope (based on effectiveA, 
where effective X is the Aapparent to a ship proceeding 
obliquely through the waves) and Z Q was faired to h at 
whatever actual X resulted (about 3.00 in the cases shown 
in Figures 2 and 3) when effective X/L-3.75 to 4.00. 

As for phase angle extrapolation, the statements 
in (2) are generally supported by vibration theory for 
a single degree of freedom system. The limiting values 
to which £ r and £_ a were extrapolated (using the same 
effective A-Hniits as described above ) were taken from 
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the following: 



£ ranges from (-180° to 0°) in head, 
fc bow and beam seas as effective A/L 
increases from zero to about 4.00. 

£ ranges from (0° to 0°) In quartering 
& seas over the effective A/L range of 
zero to 4.00. 

£ was zero throughout the 7 \/ L range for 
* following seas because of very small 
encounter frequency. 

£ ranges from (+90° to - 90°) in head and 
bow seas as effective A/L increases from 
zero to about 4.00. 

£ was 90° throughout the ?\/h range in 
” r beam and following seas. In the beam 
seas, A » ship beam and thus there 
is effectively zero pitch phase lag 
and in the following seas case u e is 
very small, 

£ ranges from (90° to 90°) in quartering 
seas • 

Before proceeding further, it is in order to pre- 
sent the basic conventions and relationships used so as 
to prevent confusion. There are many different systems 
in use so the best policy obviously was to set up a 
system and adhere to it. I have followed that policy 
based on these convent icns and relationships:* 

(a) If the wave motion amidships is 
described by 

h a h cos (^ e t ) 

the resulting pitch motion is described by 
T 0 ■ T 0 cos(^ e t +£ y ) 

*ln the notation used here a bar ( ) above a symbol im- 

plies an amplitude as distinguished from an Instantaneous 
value . 
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and the heave motion by 

Z Q sr Z D coa(o e t + £ t ) 

(b) heave is taken as positive upwards 
and 

positive pitoh angle is taken as bow 
downwards. 

(o) it can be seen from the equations in 
(a) that negative phase angles imply 
motion lagging wave and positive 
phase angles imply motion leading wave 

(d) The phase angle between heave and pitch 
is defined by the equation 

Negative values of imply heave . 

lagging pitch and positive values 
of \ imply heave leading pitch. 

(e ) Prom basic wave theory we know that 
maximum wave slope 

ZLax - ?. r k 

(f ) It can be shown that encounter 
frequency 

I O 

u>e = ^8 cos s* * 

g 

where j 2irg/^ according to basic 

wave theory. 

(g) The Dutch models were tested in waves 
of proportions: 

2 h = 1 

~ To- 

Equipped with (a) through (g) above, it is 
possible to determine the amplitude of vertical dis- 
placement at any 



*See Figure 3 for angular nomenclature. 
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( 5 ) 



point x along the ship’s length, Z x , from 
2 X 2 = 2 0 2 + X 2 Y 0 2 - 2Z Q X YqCOS li y 
where x = distance along snip’s length from amidships to the 
point under consideration, (+) forward (-) aft. In this 
thesis x = -181 feet for the 370 foot ship length and x 
changes in a geometrically similar manner for the other 
ship lengths investigated. The x distance was chosen 
to be a reasonable location of the VDS towpoint. 

The amplitude of vertical acceleration at any point 
x along the ship’s length. Ax, can be computed from 
Ax = w e 2 Zx 

The ordinate of the RAO is 

A x / g 1 2 

SrrhTXj 

and tne abscissa is In <P. 

RAO’s for ft = 180°, 135°, 90°, 45°and 0° were 
computed directly from the model data curves, extrapolated 
as necessary to cover a In w range of +0.2 to -1.6. The 
reason for the limits of this range will become evident 
late r. 

Then, in order to have RAO’s at 22.5° increments, 
the original RAO's were first plotted on an (Jp base to 
align the peaks and then crossplotted at constant values 
of £J e • The resulting RAO's f or fA = 157 l/2°, 112 l/2°, 

67 1/2° and 22 1/2 ° were then converted to a In IP base 
and plotted with the original ones, resulting in a 
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complete selection of RAO’s at 22 l/2° increments as 
shown in Figure 4. 

The RAO’s thus derived are for the full scale 
length of the destroyer hull considered, namely 370 feet 
(this was the L used in the 3./L scale of the original 
model data). 

To arrive at RAO’s for geometrically similar 
hulls of other lengths is a simple matter when they are 
plotted on a In W base (5). As stated previously, 

OJ = 1/2 

At 7)/L = a constant value, 

= k i/l/L]_ Sub script 1 indicates 

original value 

Wg = k ifl/Lg Subscript 2 indicates 

changed value 



and 




Taking the natural base logarithm of both sides 
of this equation. In - In Wg = l/2(ln Lg - In Li) and 
therefore. In w 2 - In = (f'ln CJ = -1/2 (In Ig - In I*). 
It follows that: cTln u = -l/2 lnj^j 



This last equation shows that RAO’s for other 
ship lengths can be obtained by translating the original 
RAO’s along the abscissa. 
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For example, let (fin lo = 0.1 
0.1 = -1/2 In L2 



Lo 

= .819 

L 1 

and Lg = .819 (370) = 302 » 

Similarly, if /in in = -0.1, Lg 

j— = 1.222 and Lg 



452 « 



Thus, if the RAO's are shifted to the right, keeping the 
In to scale fixed in 3pace, ship length increases and 
conversely a leftward shift causes a decrease in ship 
length. Shifting increments of 0.2, 0.1, -0.1, and 
-0.2 were chosen, resulting in RAO’s for snip lengths 
of 248 ft., 302 ft., 452 ft. and 552 ft. 

Figure 5 shows the 157 l/2° RAO shifted as 
described by the above procedure. All the others were 
shifted in exactly the same manner. 



Selection of Sea Spectra 

Through the principle of superposition, the 
large number of component waves comprising an irregular 
long-crested seaway may be represented in a spectral 
form having an ordinate representing energy and an 
abscissa representing wavelength or frequency. This is 
a well-established concept and is described rather well 
in (6). (The term long-crested is used by oceanographers 
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to describe the situation where all the waves in a system 
are coming toward the observer from the same direction 
such that he sees a well-defined pattern of long, parallel 
crests ) . 

The log-slope form of sea spectrum representa- 
tion developed by Lewis (5) was used in this present 
work. For my purpose it has one large advantage over 
other representations which have, for example, (wave 
amplitude as ordinate and encounter frequency or 
lag number as abscissa, namely; 

RAO’s for geometrically similar hulls 
of different lengths can be easily 
developed by the simple expedient of 
translating them along the In w axis. 

The ordinate of the log-slope sea spectrum is 
* *1 2 

r (lnu ) and the abscissa is lnw, 

[T/2 tT ] 

Since a statistical sampling approach was de- 
sired for this analysis, the Pierson-Moskowitz collec- 
tion of sea spectra (7), (8) was chosen because it is 
the largest, most recent collection of actually measured 
sea spectra. Pierson and Moskowitz use "corr .f t^ n * as 
their spectral ordinate and lag number, H, for their 
abscissa. See Figure 5A for an illustration of a typical 
Pierson-Moskowitz sea spectrum. The pertinent conversion 

#Corr.ft^ means "corrected” ft^-Pierson and Moskowitz cor- 
rected some spectral ordinates to smooth out occasional 
unreasonable values due mainly to vagaries in the measured 
waveheight data. 



16 



MU • It/ WU 
mmm • i * 



imi 



i*tm m >iw 



Ht. Mill 



• .911 
I .«» 

• .Ml 

• Ut 

• .Mt 
I .IN 



• MV 

• •VI 
.•V* 



• Itv 

• 111 
. 1/1 



.IN 

• IM 

• IN 



11 .ut 

u .III 
it .til 
14 .IM 
II .IN 
I* .IM 
U .IM 
M .III 
N .11/ 
*• .111 
«l . 11 * 
U .111 
II .11* 

M .IM 
*1 .IM 
M . 1*4 
♦ V .Ml 
a . i*v 
*i .in 

M .IV# 
II .Ml 
*1 •!•* 




! i 
(. ! 
1 •. I 

5 - I 

• ! ! 

i.: 

i i 
s I j 

t **| 

! I 

M 

j i 



< : 

i-l 



f !?*; 

M: i 

J N V I 

i ! 
i n ! 

i ;] 

) r- i - 



»* J 



r 



• •• «##«*» 



• i * • • • i * • • * • i • i 1 1 1 1 1 1 1 1 n i • 1 1 1 m 1 1 1 1 1 * m i • * 

>s:s:!£!:2s«:::t»:«u:x2imtu;nm}ni:»ms»s 



U)l|# IVKIU (Mwlft M#CM | Ml 



Mil • 1*1 WM 
Mt « • II 
VMM •» >M1 



.•II 

.•If 

.#11 



II 



.IM 

.III 

.Ilf 

.111 



11 .IM 
M .111 
M .IM 
M .IM 
If .IM 
M .IM 
1 * . 1*1 
M . I«f 
M .tfl 
M .IM 
11 .1*1 



4 *» 

.lit 
.Ilf 
.111 
• IM 
.111 
.IM 
. M 4 
.IM 
.IM 
.Ml 
.l*f 
.Ifl 
.IM 
.Ml 
•IM 



• 111 
.Ilf 
.Ml 
.IM 
. 11 * 




! 

f'i 

s 

i . 






s 

: *£ 

1 I 

i 

? 

-j 



0 . i tttzui :: 









/ c k£*>ejjr^T/0fiJ 

(7) t (8) 

/7 






? 

*1 

J 



i t i 



relationships are as follows: 
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H = 180 x frequency = 180 ^ 

2 TT 

g r 

= x 180 x (corr.ft 2 ) 

w ~ 

A short computer program, written in Fortran II 
source language, and shown in Figure 6, was used to con- 
vert the . Pierson-Moskowitz Spectra to the log-slope form 
used in this thesis. A sample of the print-out is 
shown in Figure 7. 

The next matter of interest is that of short- 
crestedness. This is not an easy matter to dispose of, 
simply because a short-crested sea may arise from 
different causes. As its name implies, a short-crested 
sea is one where the waves are coming toward the obser- 
ver from various directions, thus producing a short- 
crested, choppy effect. The difficulty arises from the 
fact that a short-crested sea can be caused by one 
storm system, two storm systems or any number of storm 
systems. If we place ourselves on a ship lying to in 
the ocean, maintaining a steady true compass heading, 
it is perfectly possible to conceive of a storm system 
a few miles off the bow sending wind waves toward us 
and another storm system a long distance off the star- 
board quarter sending swell waves toward us. But there 
could just as well have been just one storm system or 
more than two and they could have been off in many com- 
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binations of distance and direction, so it is only 
possible to express what is happening at a certain 
place in the ocean at a certain time by a forecasting 
or hindcasting technique. Probably the most typical 
situation is a combination of wind waves caused by a 
local storm and swell coming from a distant storm or 
storms. This is a matter of great current interest to 
oceanographers and it is my hope that in another five 
to ten years they will have accumulated a large enough 
body of directional sea spectra so as to make a rea- 
sonably sized statistical sample possible. 

However, for the moment, we must content our- 
selves with something less in our prediction methods, 
namely, either the use of long-crested seas or "uni- 
storm" short-crested seas. It is the latter that I 
have chosen. This is a situation where one storm dis- 
turbance, because of its movement anc irregular winds 
causes some spread in the wind direction and thus leads 
to a certain type of short-cre stedne ss in the seas. 

These seas are characterized by a large amount of energy 
coming from a certain direction and progressively lesser 
amounts of energy as we look away from that direction 
until there is essentially zero energy coming from about 
70 degrees on either side of the original direction. 
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This directional distribution of seaway 
energy may be characterized by multiplying a long- 
crested spectrum by a "spreading function." The one 
I have chosen is based on the work of Pierson, Neumann 
and James in (6) and is 

2 cos 2 //w 

7 r 

The result of this directional spreading of 
seaway energy can be seen at the top of Figure 9. I 
have chosen this technique because it is computationally 
practical and is at the very least, one step closer to 
reality than assuming long-crested seas. It is 
pertinent to note that my experience while working 
at the David Taylor Model Basin showed that short- 
crested predictions by this method are consistently 8- 
9 % lower than long-crested predictions of ship motion; 
at least for the head seas case. 

With the foregoing background in mind it is in 
order now to proceed with a discussion of considerations 
that led to a final choice of sea spectra to represent 
the sea severity range of interest in these predictions. 

It seems reasonable to state that as far as 
ship response is concerned, there are two characteris- 
tics of a seaway that are particularly important, namely 
the total energy, which is proportional to the area 
under the spectrum representing that seaway, and the 
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distribution of energy among the various component 
waves, an indication of which is given by the shape 
of the sea spectrum. One need only leaf through ( 7 ) 
to see that nature provides us with sea spectra of 
many different sizes and shapes. 

Now, if we wish to investigate the effect of 
sea severity on ship response, it is necessary to 
establish a primary weather severity parameter. But 
in doing so, caution must be exercised in order that 
other parameters affecting ship-response are not 
ignored . 

I have chosen significant wave height (HWj) 
as the primary sea severity parameter because it is 
proportional to sea spectrum area and thus total sea- 
way energy. (81/3 is, by definition, the average 
height of the one -third highest waves in a seaway 
record. (6)J. But, in order not to neglect the effect 
of sea spectrum shape on ship response, a family of 
sea spectra was needed where at each increment of the 
primary sea severity parameter, H^yg, there would be 
several spectra of close to the same H^yg, but having 
typical vaiations in shape among them. 

The Moskowitz family of mean synoptic sea 
spectra found in (8) is the best family presently 
available most closely satisfying the above requirements. 
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The family consists of five mean sea speotra, each 
spectrum being the mean of 8 to 14 component spectra . 
as shown in Table II. See Figure 8 for a plot of these 
five mean sea spectra. 



TABLE II 



Moskowitz Synoptic^ 
Mean Spectrum Name 


Number of Component Average of Component 
Spectra Spectra H1/3 in. ft. 


20KT 


12 


7.4 


25 


8 


13.1 


30 


12 


20.4 


35 


8 


23.8 


40 


14 


32.2 



In each subset, the H^/3 of the component spectra 
vary in most cases by no more than a few percent, thus 
assuring that the primary differences are those of 
spectral shape. 

By comparing the ship response caused by the 
H-jyj = 7.4 ft. mean spectrum with the response caused by, 
say the H1/3 = 20.4 ft. mean spectrum, one can determine 
the effect of total seaway energy on response. On the 
other hand, by comparing the ship responses due to each 
of the 12 component spectra in the = 7.4 ft. sub- 

set, for example, one is able to evaluate the effect of 
spectral shape on ship response i.e., determine the 



Moskowitz classifies his spectra on the basis of wind 
speed in knots. 
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variance in response for a particular S^y^. To my 
knowledge, it has not been attempted in previous work 
to separate these two effects and it is my hope that 
such an analysis will shed some light in our quest to 
learn more about the nature of ship response to a sea- 
way. 

The Moskowitz Synoptic Family of Mean Sea 

\ 

\ 

Spectra converted to the log-slope form is shown in 
Figure 8. It is noted that Moskowitz required that 
the sea records chosen for his family satisfy certain 
synoptic criteria as regards non-transient wind condi- 
tions, duration, fetch, etc, all of which tend to insure 
that the seas represented are fully developed or close 
to it. He has also eliminated swell from all records 
used. So one must keep in mind that the ship response 
predictions correspond to the nearly fully developed 
wind sea with no swell situation which, of course, is 
not a precise duplication of nature. 

When there is available a larger source of 3ea 
spectrum data, I believe the soundest approach would 
be to pick out all those records having very close to the 
same H^y^ without requiring that any synoptic criteria 
be satisfied. Thus, if the sample were large enough, 
one would be reasonably assured that a representative 
variation in degree of sea development would be 
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included. Each sample of records having the same 
then, would have the variation in spectral shape pre- 
viously referred to, and then one could proceed as in 
this analysis. 

Computation of Response Data 

Once the prerequisite information, namely sea 
spectra and RAO's are in hand, one may proceed to cal- 
culate ship response. 

A verbal description of the computations will 
be given to accompany the graphic illustration given in 
Figure 9. First, the long-crested sea spectrum is 
"spread" or directionalized using the 2 cos^/^w spreading 

rr 

function previously described. This operation results 
in the short-crested sea spectra shown by dotted lines. 

Next, the directional RAO's are plotted on the 
same ln^> abscissa and a decision is made as to which RAO's 
are going to be multiplied by which sea spectrum component 
at each increment of Into. Figures 10 and 10a will aid 
in this "matching" process. Say, for example, we want 
to compute acceleration response for the head seas case 
{/* = 180°) as is shown in Figure 9. We align the heavy 
wind arrow on Figure 10 with the Head Seas marker on 
Figure 10a and see that the 180° RAO must be multiplied 
by the 0° sea spectrum component, the 157 l/2 / 2Q2 l/2° 
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RAO's (v.hich are the same, of course, by symmetry) 
must be multiplied by the 22 1/2° sea spectrum com- 
ponent and so on for as many of the RAO components as 
are affected by the various sea spectrum components. 

This is done at each increment of lnu and re- 
sults in the nine response spectrum components (only 
five curves show because of symmetry). 

Now if one imagines a third dimension to the 
Response Spectrum Component curves whose axis is normal 
to the plane of the paper anc whose coordinate is ^ , 
it is possible to visualize the integration with respect 
to A of the points from the Response Spectrum Component 
curves on the crossplot for each In u> increment to arrive 
at a point on the Integrated Response Spectrum for that 
value of In Co . After the Integrated Response Spectrum 
is thus computed, one may integrate over lncjto arrive 
at the area under it which is the value of (\Te) f or the 
case considered, namely head seas, ship length = 370 ft. 
(Me is a ship response parameter that is discussed in de- 
tail in Appendix A. Briefly, the 1)E parameter reflects 
certain statistical properties of a short-term sample 
record of response during invariant conditions of weather, 
speed and heading.) The angular integration and then the 
In integration were done using Simpson's Rule in the 
computational scheme employed. Nineteen In to increments 
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of 0.1 each were used along with the 22 l/2° angular 
increments as shown. The right and left hand extremities 
of the lnvx>scale, namely, -1.6 and +0.2, respectively, 
were determined by the termination of most typical sea 
spectra on the right hand end and the termination of the 
RAO's on the left hand end. 

The process described above can then be done for 
the eight other headings between head and following seas 
so that each sea spectrum and each set of RAO’s for a 
given ship len th results in nine IfE values, one for 
each heading from head to following seas at 22 l/2° 
increments . 

This was done for the five Moskowitz Synoptic 
Lean Spectra for each of the five ship lengths and at 
nine headings, resulting in the 225 IjE values found in 
Figure 11. 

Since each one of these">fE values is the result 
of the process of computation just de scribed, where sea 
spectrum components are multiplied by RAO components to 
arrive at the Response Spectrum Components which are in- 
tegrated to arrive at the Integrated Response Spectrum 
which in turn is integrated to arrive at the ITe value, 
it can readily be seen that a very large number of 
numerical calculations are involved. However, since 
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FIGURE 11 



SAMPLE MEAN nT = ’ ACCELERATION VALLES (x) in g's 



^l/3 = ^ (Mos k°wi t z 20KT Synoptic Mean Spectrum) 



Ship 


Length — 


248 * 


302 ' 


370' 


452 * 


552' 


// 

1 


180 


.2108 


.1676 


.1238 


.0885 


.0594 


157 1/2 


.2021 


.1631 


.1212 


.0882 


.0594 


135 


.1783 


.1476 


.1148 


.0374 


.0625 


112 1/2 


.1467 


.1261 


.0999 


.0776 


.0559 


90 


.1037 


.0946 


.0812 


.0676 


.0516 


67 1/2 


.0660 


.0621 


.0544 


.0453 


.0348 


45 


.0352 


.0354 


.0335 


.0296 


.0233 


22 1/2 


.0157 


.0158 


.0151 


.0130 


.0101 


0 


.0090 


.0090 


.0090 


.0085 


.0063 



**1/3 = 


13. l' (Moskowitz 25KT Synoptic 


Mean Spectrum) 




130° 


.3181 


.2843 


.2467 


.2062 


• 


1623 


157 1/2 


.3033 


.2727 


.2369 


.1993 


• 


1574 


135 


.2644 


.24 08 


.2122 


.1815 


• 


1462 


112 1/2 


.2117 


.1967 


.1747 


.1521 


• 


1242 


90 


.1470 


.1412 


.1298 


.1168 


• 


0985 


67 1/2 


.0921 


.0901 


.0338 


.0758 


• 


0646 


45 


.0489 


.0497 


.0484 


.0452 


• 


0393 


22 1/2 


.0222 


.0224 


.0218 


.0199 


• 


0173 


0 


.0128 


.0128 


.0128 


.0124 


• 


0105 



H l/3 = 


20.4* (Mo3ko\vitz 30KT Synoptic 


Mean Spectrum) 




180° 


.3753 


.3447 


.3106 


.2742 


• 


2356 


157 1/2 


.3576 


.3300 


.2975 


.2634 


• 


2264 


135 


.3112 


.2901 


.2641 


.2363 


• 


2047 


112 1/2 


.2481 


.2349 


.2147 


.1942 


• 


1689 


90 


.1722 


.1675 


.1572 


.1453 


• 


1283 


67 1/2 


.1079 


.1066 


.1008 


.0934 


• 


0828 


45 


.0578 


.0588 


.0576 


.0547 


• 


04 89 


22 1/2 


.0265 


.0267 


.0262 


.0244 


• 


0217 


0 


.0154 


.0154 


.0154 


.0150 


• 


0131 
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FIGURE 11 (Cont'd) 



SAMPLE KEA N if E ACCELERATION VALUES (x) in g*a 

H l/S = 23 * 81 (Moskowitz 35KT Synoptic Mean Spectrum) 



Ship 



Length 

180 


— ► 248* 


302 » 


370* 


452 ' 


552 » 


.4287 


.3994 


.3646 


.3252 


.2823 


157 1/2 


.4082 


,3819 


.5488 


.3121 


.2711 


135 


.3548 


.3348 


.3085 


.2787 


.2439 


112 1/2 


.2817 


.2698 


.2496 


.2282 


.2006 


90 


.1951 


.1914 


.1812 


.1691 


.1509 


67 1/2 


.1222 


.1214 


.1158 


.1084 


.0972 


45 


.0655 


.0667 


.0657 


.0628 


.0569 


22 1/2 


.0302 


.0304 


.029 9 


.0280 


.0254 


0 


.0176 


.0176 


.0176 


.0171 


.0152 



^1/3 


= 32.2* ({•, 


loskowitz 40 


KT Synoptic 


Mean Spectrum) 


180° 


.4782 


.4524 


,4220 


.3871 


.3485 


157 1/2 


.4554 


.4323 


.4032 


.3707 


.3337 


135 


.3957 


.3784 


.3554 


.3289 


.2973 


112 1/2 


.3139 


.3038 


.2858 


.2668 


.2416 


90 


.2179 


.2153 


.2064 


.1955 


.1786 


67 1/2 


.1369 


.1366 


.1315 


.1248 


.1142 


45 


.0741 


.0754 


.0745 


.0718 


.0660 


22 1/2 


.0345 


.0347 


.0341 


.0323 


. C8 96 


0 


.0201 


.0201 


.0201 


.0196 


.0177 
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these calculations are of a repetitive nature, including 
several cyclic permutations, they are especially amenable 
to calculation by electronic computer, and this is how 
they were carried out. 

The program shown in Figure 12, based on a 
similar one in use at Webb Institute, was written by the 
author in the Fortran II source language while on temporary 
duty at the David Taylor Model Basin and the calculations 
were carried out on the UNIVAC IARC facility there. A 
sample of the print-out cm be seen in Figure 13. It 
might interest the reader to know that the computer took 
less than three minutes to compute the 225 1$E values and 
having done similar calculations with a desk computer, 

I estimate it would take about a month of 8-hour days 
to calculate the same number of “tfE values manually. 

Thus it is seen how the electronic computer cm greatly 
enhance the scope of an investigation like this one. 

The tabulation of ijE values in Figure 11 shows 
the effect of seaway total energy on acceleration res- 
ponse, but it does not show the effect of sea spectrum 
shape • 

As previously described, the effect of sea spectrum 
shape on ship response was to be evaluated by computing 
P values for each of the component sea spectra in the 
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DO 6i iO • ]»9,i 
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SHIP l 

rm$ RespoNs?~'fwr mj~s jrD“OF—; r.o is J7?ysbX^}2 

RMS RESPONSE p»R MJ SJR 0 OF 2 2 *5 TS 0* 1 «>3*E-*01 

RMS RESPONSE F.aR MJ SJR 0 OF" ’43*0 TS 0.370**F-01 

RMS RESPONSE -F*R MJ $J* 0 *F 57*5 TS 0.593«E-0t 

RMS RESPONSE Far MJ SJ* 0 m’F ~ 7 0 # 0 T S O.IOBVe 00 

RMS RESPONSE Far MJ SJR 0 0 F 112*3 TS 0.1521E 00 

RMS RESPONSE F*R MJ SJR 0 . 0? 1 3 $V3 t$ 0*1S2"E. 00 

RMS RESPONSE Far MJ • $JR- 0 OF 15 7.5 TS 0.205^>E 00 

RMS RESPONSE FaR'MJ SJR 6 mP *130.0 IS 0. ? 1 R 1 E 00 



'SHIP CeLNPfw ' 2 "JO/C. 

RMS RESPOND FaR MJ SJ3 0 OF ' 3.0 IS J.7350E-J2 

RMS RESPONSE" F^R MJ SJ^ 0 aF T?«5 !S O'.ISXaE-OT 

RMS RESPONSE F*R MJ SJM 0 OP ’4 5 • 0 IS 0.36StE-0! 

RMS RFSPO'Nsr >*RTftj £J* 0 0*“ *$> .5 TS~ 0.537TE-0! 

RMS' RESPONSE FaR MJ SJR 0 OF 70*0 TS* 0.962aE-01 

RMS RESPONSE F*R MJ SJR 6 OF H2.5 >S 0.12‘S^E 00 

RMS RESPONSr f*R MJ 5j« 0 OF 135*0 TS 0.146'lE 00 

rms RESPONsE^roR xj~§js~croF — rs rrrrs orixov. e~oo 

RMS RESPONSE F^R MJ SJR 0 O c 130.0 TS O.lS^E 00 



SHIP LtNdTH 3 37 6 

RMS RESPONSE FaR MJ SJ3 0 0* 0*0 I S' ~ J.73S0E-J2 
RMS RESPONSE F*R MJ SJ* 0 O ff 22.5 TS 0.157^E-01 
RMS R ESPONSP F oR M J 5 JR 0 OF 4 5*0 T S 0. 34 9^ E -01 
RMS RESPONSE F*R MJ S JR 0 OF 57.5 TS 0.559^E-01 
RMS RESPONSE F*R MJ SJR 0 * OF ?30 T S‘ 0.330AE-01 
RMS RESPONSE F*R MJ SJR 0 OF 112.5 tS 0.100*E 00 
RMS RESPONSE p*R‘ MJ SJ-R 0 0~F ~ 135.0 tS ‘ ~0. 1145E 00 
RMS RES°ONsp FaR MJ $JR 0 f) c 157.5 ?$• 0.120OE 00 
RMS rESPONse' FM MJ SJS 0~oF 130.0 TS 0. 1 2 l^E "OO 



SHIP Ui'MliTu “ 4 ~ • ~ - *' 

RMS RESPONSE F*R MJ SJR 0 OF 0*0 J S 0.3743E-J2 
RMS RESPONSff" F^R MJ SJR 0 MF 22.5 TS 0.134'^E-Ol 

RMS RESPONSE F^R MJ SJR 0 0 F 45*0 TS 0.?05^E-01 

•MS RESPONSE F*R MU SjR 0 OF 5 7.5 fS“ O'. '450^'E-OT 

RMS RESPONSE F*P MJ SJR 0 OF ?0*0 TS 0.532oE-01 

RMS RESPONSE F*R~MJ* SjR‘ 0" OF 112.5“ fS "0. T63*;E-0 1 
RMS RESPONSE F*R MJ SJR 0 OF 135*0 TS 0,352"5«01 

RMS RESPONSE F*R M J SJR 6 0F~ j 5 7 • S T S 0.R35«5E-01 

RMS RESPONSE Far MJ SJR 0 OF 130*0 JS 0*336^E-01 

SHIP LtNvaTu 5 £7C2J 
RMS RESPONSE F^R MJ SJR O OF^ ^tJ T $ ^ J.5627E-02 

RMS RESPONSE MJ SJR 0 O c 22*5 TS O.iOUE-Ol 

RMS RESPONSE FaR~MJ SJR 6 f)F "45*0 TS* 0*234^-01 

RMS RESPONSE F-R MJ rSjR 0 OF 57.5 TS 0.34l-?E-01 

RMS RESPONSE F^R MJ SJR 0 0*~ 70.0 TS 0."506^E-“0i 

RMS RESPONSE Far MJ_SJR 0- OF 112*5 ?s 0.522<E-01 
RMS RESPONSE FaR MJ SJR 0 OF 13*5*0 IS 0.'R72^E-01 

RMS RESPONSE F^R MJ SJR 0 oF 157*5 rs 0.503RE^01 

RMS RESPONs* Far M*J SJR 0 0 _ F 130*0 TS 0.49BSE-0T 
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five koskov/itz subsets 



Since there are a total of 54 



component spectra in the Moskowitz Synoptic Family of five 
mean spectra, converting each of these to the log-slope 
form, plotting them, tabulating and card-punching them 
as input to the ship response com outer program would 
have been too formidable a task to undertake in this 
thesis. Instead, the above was done only for the 12 



component spectra of the H 3./3 = 7.4 ft. mean spectrum 
and the 14 component spectra of the H ^/3 = 32.2 ft. mean 
spectrum so as to bracket the range of the primary sea 



straight line interpolation which turned out to be a 
quite reasonable procedure as will be discussed later. 



ship lengths in the 12 component sea spectra of the 



H 1/3 = ft. subset and in the 14 component sea spectra 



of the H -^/3 = 32.2 ft. subset were computed using the 
same computational scheme as was previously described and 
then for each combination of ship length and heading, the 
following important characteristics of the l)E distribu- 
tion were determined: 

(a) the sample mean of w as computed using 



severity parameter H 1 / 3 . The response characteristics 
for intermediate values of H- 1/3 were obtained using 



The UB values for all nine headings and five 



x = 1 



n 




I 



1/1 and xi = CtfS ) ± 



where n = 12 or 14 depending on which subset 
was under consideration and x = sample mean 
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(b) the variance of TIE was computed using: 

s%~ 

again where n = 12 or 14 and s^ = sample 
variance and s = sample deviation 

These computations also were done on the DTI" 13 
U V IVAC LARC using the Fortran II source program written 
for this purpose which is shown in Figure 14. This 
program is similar to the one in Figure 12 except that 
a section for computation of sample means and devia- 
tions was added onto the section that computed the Mi 
values. The print -out of results from this program 
may be seen in Figures 15 and 15a. 

At this point, one small issue was disposed of. 

For each ship length and heading, the mean Ml values in 
Figure 15, for example, should agree with the \Te values 
given in Figure 11 for the H ^/3 = 7.4 ft. mean sea 
spectrum, that is to say that the mean value of l(E 
obtained using the subset component spectra should be 
equal to the 1 $E obtained using the mean sea spectrum. 

This is mainly a crosscheck on the computational scheme 
used. A comparison of the data shows a maximum difference 
in the order of 4 percent, so the procedure of using a 
mean sea spectrum to obtain mean IfE response seems 
justified. Tliis conclusion engenders confidence in the 
Me values shown in Figure 11, because the ones for the 
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LARFL / '<SRl 

FGh T RAn T 'a ' 1 ’ PkOGRA’’ TO C.C"'FUIr SHIP RF SPO/SSt DATA P(_->-"> oAVkI.F 
f. vtANS A\D DEVIATIONS EEs to 

C RSI IS ! M r GA A TED ALSTONS!. SPECTrvOi)* SR 1^ slMPSONS MOL T I PL I £r< 

c SEA AN l/ RESPONSE SPEv_TKa a- >' l. AGAINST LN OMEuA 

C KC IS ^ESRORSg SP£C7A\,M Cji'VUkE.SI. 5 IS SEA SPlC t R JM 

C AREA Jj'.DtP Rbl LOCALS 0 * 0 3 3 3 3 < SuE ) » v. Is M-‘s OK sIVsl!. /«* Pi. i T JOE S 





DEV IS 


S A i v P l t DEVI 


A 1 I ON > Xd 


IS SAMP L •_ R'EA.N 




1 


F C R N ** T 


( 6 F 7 • 4 ) 










2 


F iRvAT 


! 1 C- r .„ ) 












FOR ‘•‘AT 


( 15 : 










7 


T OR vat 


I/17D S 


Hi P 


L L «\6 T *-i 


*13) 




b 


r ' ) R A T 


( 3 0 H K ,V S RE 


3P0:. 


. S L F C 


0 S J 0 ~ OF »K6»l*4r- IS > £ 1 2 • 4 ) 




9 


F OR AT 


< 1 r* i * 3 7 X * 


1 3H 


b L A b r' L 


CT^J * 2 0 X * 2 2 H ACCt-LLRAT lOr-t RESPoR 


Sc J 


lv 


f j R v A T 


(1 - * ^ 7 X • 


1 3'H 


StA 


CTRJ. * 2 0 X * IfcH V L L O C ! T v RESPONSE) 




l i 


l - 0 R ' M 


f 1 < I * ' 7X * 


I 3H 


3LA 


•:Tk • 2Ca*22h uIsRLACl • Ei\ f RESPONSE) 


4 


FORMAT 


(l’-il* A 2 X * 


2-H 


ACCt LEK 


ATl'.-N Or. VI A" IONS) 
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format 


( 1 Fi 1 * 4 4 X * 


2 *!H 


VELGC i 1 


F 0 1 V 1 m I i 0 A! S ) 




36 


FOR VAT 


( 3 AH S A 4 1_ F 


0 t V I A T i u N 


FOh .1 3 L 3 C 01 *F6*i*4H is *£14.6 


♦ 13 




in AND 


XRAk IS * F. X 4 


• 5 t 








^ H 


- OR A T 


( 1 r- * 4 ? X * 




D 1 ST'l a 0 


f. 1 . r O'. VI AT IONS) 






0 1 v<; .% s 


i ON ' *• l i 9 ) * 


-< S I ( 


iv 1 * X ( 


17*19*5) * R C * l 7 > * Ph ( 1 / ) ♦ s ( 1 9 ) * 4 


( 26 



1 • '5 * 4 ) • Q S ( ' 0 • » V ) 

A - • 5 2 4 

- = :• . 2 6 2 



C * * 


( 1 ) 


— • 


C, w 


( ? ? 


= A . 


S v 


( ^ ) 




S w 


( A ) 




CM 


( b ) 


= 7. 


S M 


( b 5 


r a * 


C M 


(7) 


^ ?. 


C V 


( F ) 


= 4 • 


zr 


(°) 


= 2 . 


CM ( 


J 0 ) 


- 4 . 


C V { 


1 1 \ 


— • 


SV ( 1? ) 


r *4 * 


s *•■ ( 


1 - ) 


- / * 


S M ( 1 4 ) 


= 4. 


SR ( 


15 ) 


r *0 • 


b v ( 


16 ) 


r a , 




17 ) 


= 2. 


S v I 


18 ) 


- 4 . 



S v ( 19) = I. 

OC is 11 = It’ 

DO 3 <5 L = 1*5 
DO 3 9 I = 1*19 

READ l*X(l* t *L)»X(P*I*i_)*X(3*I»L)*X(-*I*L)*X(5*I*L)*X!6*I»L) 

, READ 1 * X ( 7 *I*L)*X;b*I»L) •X(9*i*i_)*X( 10*1 *L)*a( 11*1*L)*X(12*I*L) 
READ 1 » X ( i 3 * I *L) *X( 14*1 * L ) »X( 15*1 *L) *X( 16*1 * L ) *X( I 7 * I *L) 

039 CONT iNUF 

READ 3* KASLS 
DO d9 1 = 1 ♦ K ASE S 

3 2 READ 2 *S( 1 > *S ( 2 ) *S< 3 ) *S( A ) *S< 5 ) *S< 6 ) *S( 1 ) *S( 8 ) *S( 5 , *o( 10 ) 

READ 2*S(ll)»Sil2)*S(13)*S<l4)*S!l3)*S<16)*S(17)*Sllo)*b(19) 

GO TO 13 3* 34 * 35 ) * El 
33 PRINT 9 
GO TO 37 
3 A DRINT 10 
GO TO 37 



fb (?-’/? A t'-jr CoofCtT Pp66r?AM TO CoA^UTf 
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O/T AX / t 6^5 



1 p ? { \; I ’ 1 

M 00 6 9 L - 1 ,5 
P R I N T 7 , L 
J'M." - 0.3027 
i,0 6 9 MUZ = ! ♦ 

0""' = j.C 

P M .Jz. = vlu'-I 
PMOZ - 22.5* PMUZ 
00 6 6 LO - 1*19*1 
PS! ( 1 Of - • 

LLC' = lO- *< 

Jf-NJ - "'ll 

*■' ' - " 1 • k 7 C " 

00 6 3 6 1 * ^ * 1 

= x.( J^u.,.0»L>*S(l ) -0.637*- , Coir ( A.Mu.% ) ) **2 

PH ( JM 1 ) =^c ( JMJ ) > J 
PJ'"J = J' o~l 

PJM S . - 22. ^ * Pj f '\i - 9u.O 

RSI l L 0 ) ~ k ■> I ( i_ ■ )+rH i J ‘u > 
iv: v,' = A - * I,- + 1 >m,i 
= - 1 

0 J TO (57,6 , -',67,o?) , N 

67 «.< JMU) - X i j ■ ■ : , lO • L I ' S ( .u i O . 6 3 7 i Cu ,M A '6 , I i '■ » 2 
PH ( jmj ) = K ( ■ ) - A 

PJ\'„ = O'' -1 

PJM'j = 72.5 * Hj'-* - 9C.3 

■ 7 5! ( L 0 ) -R - I ( 1.'., ) <-R*\ ( JVt- ! 

AM. 1 . 1/. = As " j a + 0 v ' 6 
J'Uj= 1 



Q*r '< 


CONT I i,UC 










S' ’." = RS 


I i L 


{' 1 


) * - 


1 U vJ ) * -b V. J 


066 


CONT I \ 


■ L~ 










: ( I , L • 


•U Z 
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'uRt L (5 




PRINT 


5 , 


P 




♦ 0 ! : * L f 


0 5^ 


CONTIS 


It- 










H P AO 3 


, H 










X X = < 












SO TO 


(bo 


* 


b 6 * 


vb j > L i 


56 


PRINT 


30 










GO T C 


I 2 








56 


PRINT 


31 










TO TO 


1? 










PRINT 


18 








12 


DO 16 


L r 




1 9 5 





POINT 7, i 
DO 16 MUZ = 1,9 
SUM = 0.0 

rsu-M = p.o 
1 3 DO 15 1 = 1,1? 

OS ( I, L, MUZ) = Q( I ,|_ »MUZ ) **2 
SUM - SUM + 0 ( I » L -MUZ ) 

TSUM = TSUM + US(l,L,i:UZ) 

15 CONTINUE 

XB = SUM/XK 
R = XK* ( ( XB )**2 ) 

VAR = ( 1 .0/ ( XK-: .0 ) )* ( TSUM - R) 
DEV = SORT F ( V AR ) 

N 1 = MUZ - 1 
Ml = 2 2 . 5 *N 1 
PRINT 36 » Ml , DEV, XB 



4Z 



A /(rOJPg /‘i CcH-r'p 



16 CONTINUE 



READ < 
X < = * 



GO TO ( 1 = 6 » x56 » 1 j6 ) * L 1 
156 PRINT ’0 
GO TO 88 
15b PRINT 31 
C- ) TO 8 8 
1.68 PRINT 58 
66 00 19 L = 1 o 
PRINT 7, L 
DO 1 9 Z = 1 , 9 
6 ( " 1 = .6 
To J*' = 

1 '5 :>0 17 I = I 3 » > b 

C. 6(1,1. » VL- 7 ) = 01 >Z ) *« < 

5 J'' = 5 UN + 0 ( I * L • 4U/ ) 

7 jv _ + G <• ( I , L • " 1 L ) 

1/ CuNTlNUr 

X"' = 5n'VX< 

R = y * * ( ( x ~ ) * * 2 ) 

V A 9 = II. / (>.<- - 1 . ) > * ( T 1 - a ) 

OhV r .> vRTr y VAR’ ) 

n : = >■ to’ - i 

N 1 = ^ i v 1 \ 1 

PR I N T 36 . \ 1 . Jr.V » > 3 
)? CON TIN 'E 
CmlL ^YSTE-i 



END 

TND '>8 T A b (. 
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SUB 
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SUB 
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IS 
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three intermediate values of were computed using the 

mean sea spectra only. 

The sample ifE deviations were linearly inter- 
polated over the H-jy^ range and are shown in Figures 
16 a-e for the various ship lengths and headings* 

To summarize, it is seen that at this point, 
enough information is available so that the effect on 
acceleration response of both sea spectrum total energy 
and sea spectrum shape can be examined. 

The effect of *a variation In sea spectrum total 
energy manifests itself in the five groups of sample 
mean values in Figure 11. The effect of a variation 
in sea spectrum shape manifests itself in the sample 
’{e deviations shown in Figures 16 a-e; these deviations 
from the mean exist sole ly because of this shape effect 
on ship response. 

At this point, the "\!e values in Figure 11 could 
be variously plotted to show the trends of 1$E stern 
acceleration with heading, length and sea severity, but 
before this Is done. It Is desirable to treat these T$E 
values which are characteristic of short-term response, 
as a statistical sample and use them to achieve a rea- 
sonable estimate of what $E values might be expected within 
the total universe or population of all possible destroyer 
stern acceleration "V® values during a destroyer*s lifetime 
under similar conditions of heading, length and sea severity. 
The next section will be devoted to this estimation. 
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STATISTICAL ESTIMATION OP POPULATION -tfE VALUES 



As discussed in the previous section, at each 
combination of heading, ship length and H-jyg, there is a 
mean value of ]jE (Figure 11) and a distribution of IfS 
values above and below that mean, which distribution is 
characterized by the sample deviation found in Figures 
16 a-e for the same combination of heading, length and 

H l/3* 

The assumption was made that this distribution 
is Gaussian or normal. The normal probability distribu- 
tion is completely defined by a mean value and the de- 
viation from that mean. The primary defense of this 
assumption of normality is that it does not seem un- 
reasonable in the light of present knowledge. Perhaps 
when much more work in this area has been done, a slightly 
better distribution can be found. The size of the statis- 
tical sample used here is too small to permit construction 
of a valid histogram. This assumption of normality has 
also been used by other investigators (12). 

The next and last hurdle is how to use a sample 
mean and deviation to predict the true total population 
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mean ifE and its deviation. This is where one may make 
good use of the concept of statistical inference, i.e., 
one may infer from the characteristics of a statistical 
sample the characteristics of the total population. In 
this case the assumption was made that in the sample, 
the \jE values are normally distributed with a mean,?, 
and a deviation, s. 

Sample: N (x,s^) 

In the population, the qjE values are also assumed 
normally distributed, with a meaiy^.and a deviation, (f. 

Population: N (A,<T^) 

The task now is to estimate and (J~ on the 
basis of X and s. Presently available theoretical 
statistics gives one two choices that are quite differ- 
ent in concept, namely the processes of point esti- 
mation and interval estimation. 

I will first discuss point estimation. Maximum 
likelihood estimators may be derived which give a con- 
sistent , unbiased estimate of the population parameters. 
In the case under consideration, the unbiased maximum 
likelihood estimator forAis: 

A = 1 £_ x A = x (9) 

n 

and the unbiased maximum likelihood estimator for <T is: 




**It can be shown that this expression is equivalent to 
the one used on page 40. 
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This means that the be st estimate^ one can make 
in this case, of A and <f is to use x and s, respectively. 
However this somewhat vague approach is not very attrac- 
tive unless the sample is of a respectable size, which 
in this investigation it is not. Perhaps the largest 
drawback to the use of point estimation is that when 
sample size is relatively small, one is provided with 
no specifio expression of confidence that can be had in 
the estimation of the true parameters H and (f . 

It is due to the above that I decided to use the 
concept of interval estimation. With a specific level 
of confidence in mind, say 90$, one may predict a para- 
meter interval or region which, withBI 90$ confidence, 

will cover the true parameter value or values being 

/ 

s ought • 



It has been shown (9) that the statistic 

x -A) - 

d* d/t/FT 

is normally distributed with zero mean and unit variance, 
and that the statistic 

ns 2 - 

To* 

is chi-square distributed, where do and A) are the true 
population parameter values. 

Thus, using any fconvenient set of probability 
tables, the statement may be made that 



* These tvo statistics are stochastically independent. 
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Pr f-1.15 / x -/*o / 1.15 = .875 



from the standardized normal distribution and that 



, from the distribution with n -1 degrees of freedom 
where n = 12 in this case, (n = 12 was used for the sake 
of uniformity even though the H 1/3 = 32,2 subset had 14 
component spectra-this will yield slightly conservative 
results ) 



c f */< parameter space as shown in Figure 17. 

There is a 78,8% chance that the crosshatched 
parameter region covers the true paraHaeter point (A>,(To), 
wherever it may be, that is to say that if 100 samples of 
12 each were chosen and if the 100 resulting 78.8% con- 



fidence regions were plotted as above, about 79 of them 
would cover the true parameter point and the remaining 
21 would fail to do so. This is the quintessence of 
the concept of interval estimaticn. 

The approximately 80% confidence level was chosen 

i 

as being reasonable and the following discussion is 




The joint probability is 

5^ - y * c " / — 2 ^19.*? 



(9) 



P 



£ -1.15 ^ 





(.875) (900) =* .788 



These inequalities define a confidence region in the 
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Intended to show that it is more than conservative 
enough for the purpose of these estimations. 

The maximum likelihood estimators tell us that 
P is the most likely location of (/b, (To). With this 
and the approximately 80% confidence level in mind, I 
have chosen (Ai, <5^u) which will give the ’'worst” normal 
distribution in terms of extreme values as the true 
parameter point. This is a fairly conservative choice, 
beoause of the likelihood that (<fo,A>) is somewhere 
near P. I feel this is a situation where a little 
judicious conservatism is better than being inordinately 
bold in ignorance . 

The parameters of N(Ai, 6 u^ ) may thus be evalua- 
ted from the following simple relationships: 

/*u s* x + .540 s 
<fu = 1.620 s 

These equations show the amount by which the 
sample parameters are increased so as to arrive at an 
estimate of the population parameters. 
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INVESTIGATION OP TRENDS IN ]Je ACCELERATION 



From this point, mean ifE acceleration will imply 
estimated population., mean ~l|E acceleration ( /*u) and 
acceleration deviation will imply estimated population 
deviation of / \|E acceleration (<£u). 

From the data in Figure 11 (x) and Figures 16 
a-e (s), many combinations of ship heading, length and 
H ^/3 can be used to investigate trends in accelera- 
tion. For purposes of the present work, I have considered: 

(a) one ship length (570 ft) - all headings- 

all § 1/3 

(b) one heading (head seas) - all lengths- 

all H l/3 

These cases show the effect of ship length and heading 

x/' — _ 

as well as on mean uB acceleration response and its 

deviation. Other cases could, of course, be considered 
if desired, using other data from Figures 11 and 16 a-e. 

Figures 18 and 18a deal with the 370-foot ship 
length case and show how mean ifE acceleration and its 
deviation vary with heading and H-jy 3 . Figures 19 and 19a 
show trends with ship length and H 1/3 for the head seas 
c & s© • 
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Finally, by crossplotting Figures 18 and 18a 
at discrete values of H 1/3 (which can be thought of as 
different sea states - see Appendix B), Figure 20 was 
derived. Similarly crossplotting Figures 19 and 19a led 
to Figure 21. In the interest of clarity, only enough 
deviations are plotted on Figure 20 to demonstrate the 
fact that the deviations hehave in much the same fashion 
as do the mean ^ values for any given ship length, i.e., 
a plot of deviations results in curves similar to the 
curves of mean acceleration in Figure 20 . 

If one accepts this conclusion of similar be- 
havior of means and deviations, the fairly mild curva- 
ture of the mean iJE curves between Hjyg = 7.4 and 32.2 
feet in Figure 18 tends to justify the assumption of 
linearity used in interpolating the sample deviations in 
Figures 16 a-e in the same H-jy^ range. 

There are some important points to be learned 
from an examination of Figures 20 and 21 . 

Figure 20 illustrates a point of great opera- 
tional significance o It is seen that if destroyers were 
designed so as to tolerate beam, quartering and following 
seas operation from a rolling and maneuverability stand- 
point, the stern accelerations could be greatly diminished 
and thus the effectiveness of the DD/VDS system is 
greatly enhanced in any given sea state. A step in the 
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right direction can be taken in the installation of 
anti-rolling fins on destroyer-type ships, because as a 
consequence of being able to tolerate beam seas opera- 
tion, Figure 20 shows they could diminish stern accele- 
rations to between 50 and 60 percent of the head seas 
accelerations . 

As a slight digression, I see no reason why 
the same reasoning cannot be applied to problems associa- 
ted with bow movement such as slamming, shipping of green 
water, degradation of bow sonar performance and general 
ability to maintain high speed in rough seas. 

Figure 21 illustrates the fact that ship length 
has a greater percentage effect on stern acceleration in 
small waves than in large ones. Consequently, increased 

I 

ship length pays significant dividends during the majority 
of the time because the small waves are much more prevalent 
than the large ones, but pays lesser dividends when con- 
sidering relative performance in severe waves. 



* 
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UTILIZATION OF tJE ACCELERATION DATA* 



A crucial question is how to most effectively 
make use of data such as that plotted in Figures 20 and 
21, By the means discussed in Appendix A, it is possible 
to predict from the l)E value the probability of occurrence 
of any particular value of stern acceleration. 

However* it is highly desirable to simultaneously 
consider the normal distribution of "tfE values (and thus 
the Rayleigh parameter R) and the Rayleigh distribution 
of acceleration single amplitudes. The theory for this 
has been developed by Nordenstrom in (12), He assumed 
that the parameter R in the short term Rayleigh dis- 
tribution is normally distributed and proceeds to obtain 
overall distributions by numerical integration of the 
convolution of the Rayleigh and normal distributions. A 
sample of the results of this numerical integration can 
be found in Figure 22 with suitable nomenclature changes 
made to adapt it to the present work. One can enter this 
diagram with (ju/Z^u and for different values of N (number 



‘fit . 

A detailed discussion of the ijE parameter is given in 
Appendix A. 
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of cycles) obtain a factor, which when multiplied by^u 
allows one to derive curves of Ax/g such as those shown 
in Figure 23 for the head seas and 370 ft ship length 
case . 

It must be pointed out here that although these 
curves in Figure 23 could be termed short-term distribu- 
tion curves in accordance with Jasper f s definition (11), 

I will not do so, because it is absurd to consider that 
millions of cycles could take place during a short term 
of say, twenty minutes. However, I do not choose to call 
them long-term distribution curves, either, because this 
implies a natural vaiation of ship speed and heading, 
which is not the case here. These curves will be con- 
sidered "long-term” only in the sense that each curve 
is the sum of many short-term distributions, keeping in 
mind the assumptions of invariant ship speed and heading 
that I have been forced to make, the former because of a 
paucity of model test data and the latter because of a 
lack of information available as to probability distribu- 
tion of heading to the sea. This matter is discussed 
further in the next section. 

A partial answer to this dilemma is to consider 
the probability of occurrence of the five sea states 
shown. It is then possible to weight these five curves 
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by the probability of occurrence of each sea state and 
arrive at a single curve that is at least one step 
closer to being a true long-term distribution curve » 

The following sea state distribution from (14) 
for the North Atlantic (40-50° North Latitude ) was used 
as an example : 

H^/3 in feet Probability of Occurrence 



Calm 


1.3% 


less than 4.92 


32.7 


4.92 - 9.84 


44.1 


9.84 - 14.77 


14.9 


14.77 - 19.69 


3.9 


19.69 - 26.20 


2.1 


greater than 26.20 


1.0 

Too7o% 



Plotting this information on probability paper 

>A 

gives one a continuous distribution of H1/3. By breaking 
this continuous distribution up into discrete increments 
centered on the H1/3 values being used for prediction 
purposes (6,11,17,25 and 37 feet), it is reasonable, but 
not mathematically rigorous to state that: 

&l/3=6 ft is representative of what occurs 70.0% of the time 



n 


11 


91 


22.0 


II 


II 


17 


II 


5.5 


H 


99 


25 


99 


2.0 


99 


It 


37 


99 


.5 


II 








I60.0% 





These statements are not mathematically precise 
because the ship response due to a discrete ^1/3 is said 
to be equivalent to the ship response resulting from a 
range of 1^/3 values above and below that discrete value. 
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and this is obviously affected by the size of the H-jyg 
increment between discrete values of H^/ 3 * The larger 
the increment, the larger would be the deviation in 
values; larger than the deviations caused merely by 
variation in sea spectral shape. 

However, in the analysis undertaken here, the 
degree of approximation seems reasonable with general 
regard mm to accuracy of the predictions. 

The five curves in Figure 23 were weighted at 
various levels of acceleration by the above probabilities, 
and the result was the curve shown in Figure 24 which I 
choose to call a quasi-long term distribution curve for 
a particular length ship under the assumed conditions 

* r 

of ship speed and heading. When statistical information 
is available regarding speed and heading distributions, 
inclusion of this information in the analysis will make 
a curve such as in Figure 24 completely general, and 
then it could honestly beccalled a long-term distribu- 

' '4 . 

tion curve ** 

Fortunately, however, a quasi-long term distri- 
bution curve such as. that in Figure 24 is presently 
useful for systems analysis and design because the con- 
ditions of speed and heading used in this work are such 
as to yield conservative results. In other words, more 
work and increased knowledge will show us how much 
lower our predictions can safely be. Here again this 
is given more disoussion in the next section. 

i 
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A DEMONSTRATION OP THE SYSTEMS VIEWPOINT 



With predictions available such as those 
developed in the preceding seotion, it is possible to 
make an evaluation of system effectiveness, keeping in 
mind the assumptions as to ship speed and heading that 
were made. The significance of these assumptions will 
be discussed further later. 

In the case under study, the system is the 
DD/VDS system. Destroyer vertical stern acceleration 
is no doubt the most important single parameter affecting 
operation of the DD/VDS system, given an otherwise oper- 
able system. The structural ability of the cable to 
remain intact in rough weather is a go/no-go situation 

4 $ 

dependent on vertical stern acceleration. Vertical 
stern acceleration also bears on understanding and 
solving the VDS control problem and on the design of 
the towing device. 

A systems analysis has been done herein using 
vertical stern acceleration as the significant parameter- 
other parameters could be included in the same manner, 
if 8uchw*e deemed necessary. 
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There are two basic areas of interest in this 
analysis of DD/VDS operational effectiveness based on 
vertical stern accelerations; they can best be posed as 
questions : 

(a) in the case of an extant system, where a 
design assumption as to vertical stern 
acceleration has been made, what should 
the upper operational limit be, with 
respect to weather, in order to give 

« 

reasonable assurance that the cable will 
not be snapped and that the VDS might be 
controllable, that is to say, in order 
that the design assumption of vertical 
stern acceleration will not be exceeded? 

(b) in the case of a system under original 
feasibility or design study, what value 
of vertical stern acceleration should be 
used as a design requirement, so as to 
insure that the system will be opera- 
tionally effective for an acceptable per- 
centage of the time on a long-term basis? 

For the extant system, the curves in Figure 23 are 
the most useful. Say, for example, that Ax = 1.5 g’s was 
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used as a design assumption and therefore the design was 
based on this figure. Let us say, further, that it is 
desired that the probability of exceeding 1.5 g be .001% 
or less. Figure 23 tells us that in seas characterized 
by H ^/3 = 25 feet, the probability of exceeding 1.5 g is 
slightly less than .001%, which is satisfactory. However, 
in the case of seas characterized by H-jy 3 = 37 feet, the 
probability of exceeding 1.5 g is about .095% , which is 
too great. Therefore, we would specify that this system 
avoid operating in seas having H 1/3 greater than about 
25 feet. Since seas of this severity might be expected 
to occur only two or three percent of the time, this 
system could be expected to be operationally useful 
97-98 percent of the time, considering the restriction 
imposed and the probability of failure that was 
specified as being acceptable. The numbers used in this 
example are for illustration only - if a smaller prob- 
ability of failure is desired, for example, a lower 
upper sea severity limit would have to be imposed, thus 
reducing the total percentage of the time that the system 
could be operated. This is obviously a situation involving 
various trade-offs and one must pick the combination best 
suited to his needs. 

In the case of a new design, it can be seen that 
Figure 24 is the most useful, if it is desired, as it 
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certainly seems it would be, to avoid any sea severity 
upper limitation on the operation of the system. Once 
an assumption is made as regards likely areas of opera- 
tion, a curve such as that shown in Figure 24 can be 
used to arrive at a reasonable value of stern accelera- 
tion for use as a basis for design. 

For example, let us consider a DD/VDS system 
that is at sea for 40 percent of a ten-year life span 
(the ship will last longer, but most, likely the weapons 
system will be changed), and assume that the VDS is 
streamed for half the time tile ship is at sea. One 
could expect about 9 million cycles of stern accelera- 
tion during that period of time, based on en overall 
period of stern motion of about 7 seconds. Entering 
Figure 24 with N = 9 million, one rends off about 1.8 g's 
as the level of stern acceleration above which we can 
expect only one of the 9 million cycles to occur. The 
corresponding probability of failure if we use 1.8 g's 

T 

as a basis for design is about .00002 percent. 

It must be emphasized again at this point, as in 
the previous section, that curves similar to those being 
discussed are only as good as the assumptions made in 
deriving them. In the specific case under consideration 
in the present work, the curves in Figures 23 and 24 are 
based on a single ship speed of 17 knots and, of course. 
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the system will not spend its entire lifetime operating 
at 17 knots, although some weight may be given to the 
fact that 17 knots could be considered a typical speed 
for operation of the DD/VDS system. Another fact that 
must be kept in mind regarding these particular figures 
is that they are for head seas operation and, here again, 
the system will not spend its lifetime operating in head 
seas. As a matter of fact, as discussed previously, it 
would considerably lower the levels of stern accelera- 
tion if we were to design our destroyers with anti- 
rolling devices so as to be able to avoid head seas 
operation. Therefore, the curves given here will yield 
conservative results in any case. 

My primary purpose in this section is to present 
a method for systems analysis that can be quite useful 
in evaluation of system effectiveness when the input to 
the analysis is reliable in magnitude and sufficient in 
scope. At present, I consider the available input in- 
formation marginal in some respects, namely: 

(a) probability distributions of ship 

heading to the seame not presently 
available and should therefore be 
given some attention. It seems evi- 
dent to me on the basis of four years 
at sea on board destroyers that as 
weather becomes more severe, the 
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tendenoy to steer so that seas are on 
the bow increases because of a desire 
for minimum motions and maximum 

i 

maneuverability. It is due to the 
foregoing that I assumed head seas as 
being slightly more conservative than 
an assumption of bow seas. However, 
the more widespread use of such de- 
vices as anti-rolling fins will no 
doubt change the situation. 

(b) it would be highly desirable to include 
the effect of speed in these predictions 
even though it probably has a second order 
effect compared with ship length and sea 

¥ 

severity. At present there is no model 
test series available that covers all 
headings and speeds. Perhaps analytical 
derivation of the RAO's, such as is now 
being investigated at Webb, would solve 
this dilemma, by replacing the model 
testing process with an analytical pro- 
cedure than can be done with a computer. 

One further point that must be discussed is that 
of attempting a complete systems analysis. There are 
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many factors that affect a system as complex as a ship 
and its equipment, and I have examined only one sub- 
system of the larger overall system. The reliability 
of the machinery in its task of motivating the ship 
during VLS operations at the desired speed certainly 
affects the effectiveness of the DD/VhS system, but the 
inclusion of factors such as this is left to the opera- 
tions analyst charged with maintaining the macroscopic 
view, more commonly referred to as the "big picture." 

It is my hope that the information developed in this 
thesis regarding destroyer stern acceleration might aid 
him in his task. 
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CORRELATION OP PREDICTIONS WITH FULL-SCALE MEASUREMENTS 



As a result of interest in VDS systems, the U.S. 
Navy Underwater Sound Laboratory in New London, Connecti- 
cut, has recently embarked on a program of collecting a 
statistical sampling of vertical stern accelerations on 
destroyers at sea (17). For this purpose, Giannini 
Controls Corporation Model 2432 statistical accelero- 
meters have been installed on various destroyers at the 
stern, as well as other locations. 

Up to the present time, these accelerometers have 
been installed so as to detect instantaneous values of 
positive acceleration that exceed the following four 
thresholds: 1.20g, 1.35g, 1.50g and 1.65g, and record 

on the four counters, the number of exceedences above 
each of the thresholds,, 

It must be pointed out here that since these 
accelerometers are designed to be used horizontally as 
well as vertically, the zero reference for neasuring 
vertical accelerations is one g rather than zero g, the 
latter being used throughout this thesis. As a oonse- 
quence, the above thresholds become .20g, .35g, .50g 
and .65g in the zero g reference system for vertical 
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accp lerations . No attempt will be made here to debate 
the pros and cons of the two systems - suffice it to say 
that there are arguments for both. 

As of the present, there has been collected a 
small amount of this statistical data (and to my know- 
ledge it is the <?nly such data available) and it is of 
considerable interest to compare it with the predictions 
made in this thesis. 

The data received from the Underwater Sound 
Laboratory is as follows: 

USS UOALB (DD693) 

Number of accelerations above .20g = 708 

.35g = 239 
.50g = 44 

. 65g = 4 



Total time of test = 226 hours 
Estimated number of cycles for 6 second 
period of motion = 135,600 cycles 
Weather 1/6-6 foot wave heights observed 
Headings and speeds - various 

The probability of exceeding various levels of 

acceleration, therefore, would be as follows: 

Probability of exceeding .20g = 708/135,600 = .530% 

.35g = 239/135,600 = .176 

. 50g = 44/135,600 = .0325 

. 65g = 4/135,600 = .00294 



USS PUTNAM (DP 757) 

Number of accelerations over .18g = 6009 

.38g = 1529 
,4 6g = 149 

. 67g = 17 
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Total time of tost 612 hours 
Estimated number of cycles for 6 seoond 
period of motion = 367,000 cycles 
Weather 1-10 foot waveheights observed 
Headings and speeds - various 

Probability of exceeding .18g = 1.64% 

. 3 8g = .43 

• 46g = .041 

. 67g = .0046 



USS H. J. ELLISON (DP 864) 

Number of accelerations above .19g = 1229 

•34g = 156 

,48g = 18 

• 6 6g = 1 

Total time of test 629 hours 

Estimated number of cycles for 6 second 
period of motion = 377,000 cycles 
Weather 0-12 foot wave heights observed 
Headings and speeds - various 

Probability of exceeding .19g =.34% 

.34g =.043 
,48g =.0049 
. 66g =.0003 



The Moale and Putnam are destroyers of the Sumner 
Class whose length on the waterline is 369 feet. The 
Ellison is a Gearing Class Destroyer whose length on the 
waterline is 383 feet. 

The probabilities of exceeding the four levels 
of acceleration for each ship were plotted in Figure 23 (Pg. 70 ) 
to compare them with predicted values. 

Since the predicted curves are for head seas 
and 17 knots ship speed, one would expect them to be 
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high relative to the observed data which is for various 
headings and speeds. Therefore, since the observed data 
was taken in seas in the neighborhood of 0-11 feet in 
height, one would expect the observed data to plot some- 
where below the - 11 feet curve. As can be seen in 

Figure 23, this is exactly where the observed points plot. 
The dashed line shows the trend for the two 3 69 foot 
ships and the dot-dash line shows the trend for the 383 
foot ship. The longer ship exhibits a tendency toward 
lower stern accelerations which agrees with the trends 
of stern acceleration with ship length previously dis- 
cussed . 

The faot that the lines drawn through the ob- 
served data have a larger slope than the individual 
curves seems to reflect the fact that the observed data 
included the effect of a range of weather conditions. 

The curve of predicted values in Figure 24 has a larger 
slope and it likewise includes the effect of a range of 
weather conditions. 

This section is intended in no way to be a com- 
prehensive correlation of observed data with predicted 
data. Indeed, much more observed data is needed so one 
might cull it for information as to accelerations for a 
particular heading and speed. However, the information 
presented here tends to increase one's confidence in the 
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quantitative and qualitative aspects of the stern 
accelerations predicted by the method presented in this 
thesis . 



* 
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CONCLUSIONS 



1. The method presented in this thesis pro- 
vides the sys tem designer with a powerful technique — 
where ship response to the sea is a pertinent factor— 
for analyzing on paper the worth of a proposed system while 
still in the conceptual stage*' 

2* The available full-scale stern acceleration 
data tends to confirm the qualitative and quantitative 
aspects of the predictions made by the method presented 
herein* 

3* Prom the standpoint of diminishing stem 
accelerat ion, ship heading to the seas is the most impor- 
tant single controllable factor and, therefore, destroyer- 
type vessels with installed equipment sensitive to stern 
accelerations should be designed so as to be able to 
tolerate rough seas operation in beam, quartering, and 
following seas* Installation of anti-rolling devices is 
a step in the right direction by enabling the ship to 
operate in beam seas and thereby paying a dividend of 50- 
60 percent diminution in stern acceleration* 

4* Ship length has a greater percentage effect 
on stem acceleration in light seas than it does in 
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heavy seas. Consequently, Increased ship length pays 
significant dividends during the majority of the time 
since the light seas are much more prevalent than are the 
heavy ones* The effect of ship length is not so marked 
when considering relative performance on a percentage 
basis in severe seas* 

5* Sea spectrum shape has a significant effect 
on stern acceleration response and should therefore be 
accounted for in the method used to predict that response. 
The method used here seems to be effective in that regard* 

6* As demonstrated in this thesis, the system 
approach seems particularly useful in certain aspects of 
ship design, since a ship-related system may most logically 
be considered in terms of its intended mission, for example, 
the location of enemy sii>marines* 

7* In a problem such as the one studied here, 
the use of a high speed electronic computer virtually 
frees the designer from worries regarding time required 
for calculations and allows him to spend full time in- 
suring the soundness of the analytical processes being 
used* 

8* The generality of the analysis developed 
here insures that a similar method can be utilized to 
investigate other system characteristics related to sea- 
keeping, such as shipping water over the bow, sonar de- 
gradation and slamming* 
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9. The analysis used in this thesis shows that 
the behavior and effectiveness of a ship-related system 
is ideally suited for expression in probabilistic terms. 
However, both ship designers and operators will need to 
develop a whole new framework of thought and reference 
as to acceptable probability levels of system failure 
and/or loss. Indeed it puts a designer face to face 
with his responsibilities when he has to express loss of 
resources, human and otherwise, in mathematical terms 
that are uncomfortably precise. 
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RECOMMENDATIONS FOR FUTURE WORK 



1. Investigate means that would allow destroyer- 
type vessels to tolerate operation in rough quartering 
and following seas. The obstacles to be overcome include 
considerations of stability and maneuverability. 

2. Continue work presently being done in the 
development of theory for analytical derivation of res- 
ponse amplitude operators, including experimental verifi- 
cation. 

3. Investigate probability distributions of ship 
heading to the sea under various conditions. 

4. When data becomes available, undertake a 
comprehensive correlation between full scale stern 
acceleration measurements and the predictions contained 
herein. 

5. Give increased emphasis to collecting more 
seaway measurements useful for ship design, preferably 
directional data measured with a buoy rather than a 
shipborne recorder. In this regard, effective liaison 
between ship designers and oceanographers Is of the ut- 
most importance. 

6. Investigate stern vertical velocities and 
displacements under various conditions of weather, ship 
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length and heading using a method similar to the one 
presented here. Also investigate the effect of ship 
speed on stern vertical acceleration, velocity and dis- 
placement . 

7. Investigate the application, to the VES 
cable stress situation, of an instrument similar to the 
flexing stress monitor described in (18). This instru- 
ment could be used to compute and present maximum probable 
cable stress during the ensuing 15 minutes or so of 
operation and thus give the ship's personnel an objec- 
tive basis for making speed and/or heading changes when 
the probable cable stress becomes excessive. 
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APPENDIX A 



DISCUSSION OP THE lE PARAMETER 
A reoord of any wave-caused phenomenon! in 



irregular seas such as stern acceleration would look 
something like 




where xv is stern acceleration measured from the mean 
value at equal or random time intervals and y A is peak to 
peak stern acceleration measured wherever it occurs. 

It is now widely accepted in the literature* that 
the x/ ? s are normally distributed and that the y^’s are 
Rayleigh distributed on a short-term basis. A short- 
term distribution is applicable so long as uniform condi- 
tions of sea, ship speed and course persist. This is to 



*See for example, (11) and (13). 
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be distinguished from a long-term distribution which allows 

the sea, ship course and speed to vary over a range of 

conditions . Jasper introduced these definitions in (11) 

and they are used by many investigators, including Nordenstrom (12), 

(Peak to peak) Rayleighs ^ (y) • 2jr e -y^ A 

A 

(Sucessive r .(*i - //) 

Ordinates) Normal: f (x) • 1 . 2 a~* * 

J <rTW * 

It is seen, then, that the parameters necessary to define the 

2 

above distributions for a particular case are R and 0~ , Fortunately, 
it is possible to get these parameters from the Integrated Response 
Spectrum which, through the superposition concept, has certain 
mathematical properties in common with the irregular record shown. 




* for brevity, the normal probability distribution 
function is often referred to as N(^// 
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It has been shown in the literature* and is now 
widely accepted that, for the type of spectrum being used 
in this work. 



R = 8<T 2 = 4E (13) 

Thus it is seen that R and C 2 are known once E is known 
from the Integrated Response Spectrum. 

It can also be demonstrated, using the concept 
of maximum likelihood estimation, that 




and since root mean square (of peak to peaks) is defined 



as 



/peak to] 
RMS 1 peaks ) 



V / * 0 ^ 

' ' u n 

one may relate RMS( pt p) to R and thus E as follows 
RMS/ 



(pt p) ■ t Ljf . =im = 2® 



RMS of peak to peaks is commonly referred to as RMS of 
double amplitudes. Therefore RMS of single amplitudes 



RMS/single \ * RMS . /2 = ifB 
(amplitude^ pt;p 

Thus it is seen that the if E values plotted in 
the various figures can be thought of as being the root 

*5>ee (13) 
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mean square of single amplitudes of stern acceleration 
and such is the case throughout this thesis. 



All of the above relationships were originally 



developed using an irregular record of wave height and 
the corresponding sea spectrum, but I have discussed 
the situation in terms of stern acceleration because 
ship characteristics dependent upon waves such as pitch 
and heave have been found to behave in the same fashion 
as the waves. 



Since most often one is most interested in single 



or double amplitudes of motion, it is common practice 
to use the parameter B which defines the parameter R 
and thus the Rayleigh distribution of amplitudes to 
arrive at, for example, the average highest expected 
value of stern acceleration single amplitude in 20 cycles 
from 



N>50) may be computed to within 3% accuracy using the 
following formula from (10): 



Examples are as follows: 

N=number of cycles Average of highest expected 
values of single amplitude 




Multipliers for other numbers of cycles N (where 



(In N) 1//2 + x (In N)“ 1//2 
2 



where T = Eulers Constant = 
.5772 



50 

100 

200 

500 

1000 



2.12 7/B 
2.28 TfE 
2.43 -tfE 
2.61 TtE 
2.74 TfE 
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the following 



For the reader’s convenience, 
are also quoted : 



( 6 ) 



Most frequent single amplitude 
Average of all single ampli- 
tudes 

Average of l/3 highest single 

amplitudes 

Average of l/lO highest single 

amplitudes 



0.707 if E 
0.886 ifB 
1.416 ZJB 
1.800 fS 
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APPENDIX B 



A UNIFORM DEFINITION OF SEA STATES 
For about fifteen years, a world -wide sea state 
scale (World Meteorological Organization Code 75) (15) 
has been in general maritime use by both naval and 
merchant ships for the collection of weather data at 
sea. 



It seems logical, therefore, in the interest of 
meaningful full-scale correlation, that all our engineering 
predictions for ultimate application to the ships at sea 
should be in terms of WMO Code 75, which is quoted below 
for the reader's convenience: 



Code 

0 

1 

2 

5 

4 

5 

6 

7 

8 
9 



STATE OF SEA -WIND WAVES (WMO CODE 75) 



Description 

Calm-glassy 

Calm-rippled 

Smooth-wave le t s 

Slight 

Moderate 

Rough 

Very Rough 

High 

Very High 
Phenomenal 



Height in Feet * 
0 

0-1/3 
1/3-1 2/3 
1 2/3-4 
4-8 
8-13 
13-20 
20-30 
30-45 
over 45 




t is proposed herein to consider this oolumn as being 
1/3. 
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In apite of the above, other aea state scales, 
such as the one reproduced In Figure 25, are in wide use 

e 

by designers, presumably because of the mathematically 
convenient use of significant wave height (Hjyj) versus 
sea state designations. However, one simple and defenslbls 
assumption would remove the need to use anything but WMO 
Code 75, This is discussed in the two following paragraphs. 
Both references (6) and (16) agree with what many 
oceanographers seem to feel in that at-sea wave height 
observers tend to overestimate wave height. Reference 
(16) even goes so far as to instruct the observer to 
record the average of the estimated heights of the 
larger, well-formed waves. 

Consequently, it seems to make a great deal of 
sense to specify the assumption that observed wave height 
equals significant wave height, H^, where is the 

average height of the 1/3 highest waves in a seaway. This 
would seem to bridge the gap between the designers who 
have been prone to use such scales as the one shown in Figure 
25 because of its mathematically desirable use of and 

the ships at sea who have been making observations based 
on WHO Code 75, 

It might be pointed out that the values of 

6, 11, 17, 25 and 37 feet used in this thesis for pre- 
dictions of stern acceleration wepe chosen as being 
representative of WMO Code 75 Sea States 4, 5, 6, 7 and 
8 respectively. 
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